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ANALYTIC CROSSING PROBABILITIES FOR CERTAIN 
BARRIERS BY BROWNIAN MOTION 

By Nabil Kahale 
ESCP-EAP 

We calculate crossing probabilities and one-sided last exit time 
densities for a class of moving barriers on an interval [0, T] via Schwartz 
distributions. We derive crossing probabilities and first hitting time 
densities for another class of barriers on [0, T] by proving a Schwartz 
distribution version of the method of images. Analytic expressions for 
crossing probabilities and related densities are given for new explicit 
and semi-explicit barriers. 

1. Introduction. Let Wt be a standard Brownian motion and g a con- 
tinuous function on an interval [0,r]. For any real number u, consider the 
first hitting time 



(1) Tin) 



in{{t e[0,T):u + Wt= g{t)}, if such t exists, 
T, otherwise. 



Define the boundary-crossing probability with respect to g as -P(r(0) < T). 
The calculation of boundary-crossing probabilities or other functionals of 
the first hitting time density arises in several areas such as statistical test- 
ing [7, 23], the valuation of barrier options [9, 14, 24] and default mod- 
eling [2, 11]. Closed-form formulae for boundary-crossing probabilities or 
first time hitting densities are known for few barriers other than the well- 
known case of linear boundaries. For quadratic boundaries, analytic rela- 
tions [10, 25] between the first hitting time density and Airy functions 
have been established. First hitting times for square-root boundaries have 
been studied in [1, 7, 26]. Analytic expressions involving hypergeometric 
functions are given in [17] for densities of the first hitting times mf{t > 
0:Wt>a + b^/c + t} and inf{t >0:\Wt\> b^/c + t}. Using the method of 
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images [4, 5], crossing probabilities and first hitting time densities have 
been calculated for a class of moving barriers via a-finite measures. The 
method of images is studied in detail in [15], Chapter I. Using the well- 
known result that tWi/^ is a Brownian motion, a correspondence has been 
established [15], Chapter I, between the method of images and the so-called 
method of weighted likelihood functions [22, 23] (see also [12], Section 4.3). 
A relation between the densities of the last exit time and the first hitting 
time for the time-reversed boundary is shown in [25]. 

This paper provides crossing probabilities for new moving barriers. Our 
approach is inspired from the reflection principle (see, e.g., [21], page 105) 
that shows that a constant barrier-crossing probability is twice the proba- 
bility that Wt is beyond the barrier. The reflection principle can be easily 
generalized as follows. If g{0) > and A is a Borel subset of [g{T), +oo) 
such that the probability that a Brownian motion starting at g{t) at time t 
belongs to A at time T is a constant c independent of t, then P(r(0) < T) = 
P{Wt G ^)/c. Theorem 2.5 gives a further generalization of the reflection 
principle by showing that the crossing probability for certain moving barri- 
ers can be retrieved from the density of Wt via Schwartz distributions. The 
probabilistic basis for the proof of Theorem 2.5 is closely related to Fortet's 
equation ([8], page 217), 

VT \ VT J \^T-t{u) \^T- t{u) 

valid for u < g{0) and v > g{T), where N{x) = (27r)-i/2 J^^ exp(-nV2) du. 
Fortet's equation has been used to establish integral equations [18] and nu- 
merical approximations (see [2] and references therein) for the first hitting 
time density. To our knowledge, though, it has not been previously applied 
to establish closed-form formulae for curved boundary crossing probabilities. 
Section 3 gives an analytic expression for the one-sided last exit time density 
for the barriers considered in Theorem 2.5 under an additional smoothness 
condition. In Section 4, we show through time-inversion that results in Sec- 
tions 2 and 3 yield hitting probabilities and first crossing-time densities 
for certain boundaries on [T, -|-oo). Section 5 generalizes our results to two- 
sided boundaries. Section 6 establishes a Schwartz distribution version of the 
method of images by building upon our proof techniques. Finally, Section 7 
contains concluding remarks. The classical and the Schwartz distribution 
versions of the method of images allow the calculation of the crossing prob- 
ability and of the first hitting time density for a class of barriers which is, 
to our knowledge, unrelated to the one obtained via Theorem 2.5. Several 
examples illustrate our results. Except for one-sided and double-sided linear 
barriers, where crossing probabilities and related densities are known or can 
easily be derived through classical tools, the results in our examples are new. 
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to the best of our knowledge. In particular, we give simple explicit formulae 
for boundary-crossing probabilities and the one-sided last exit time densities 
when g{t) = a + b\/T — t and, under suitable conditions on a and b, when 

g{t)=a-^{T-t)H^^). 

2. Crossing probabilities. Define the function H from to by 
(2) H{u,v) = P{t{u)<TAu + Wt>v). 

If /i is a Schwartz distribution or a (T-finite measure and u G M, let 

u — V 



U{u,t;fi) 




dfi{v), ifO<t<r, 
otherwise. 



A technical difficulty often encountered in this paper is the interchanging of 
Schwartz distributions and integration in the measure-theoretic sense, which 
is valid if certain conditions are met, as shown in Lemma 2.1, but may lead 
to erroneous conclusions if they are not, as shown in Example 3. 

Lemma 2.1. Let mi be a finite measure on m2 a Schwartz distribu- 
tion on M with compact support K, and (^{t, v) a measurable function on R^. 
Assume there exists an open set V containing K and a sequence such 
that, for any A; > 0, -^4){t,v) exists and is bounded in absolute value by 
for {t,v) £RxV. Then: 

1. The function v J (j){t,v) dmi{t) is C°° on V. 

2. The function t^ J (p{t,v) dm2{v) is bounded and measurable on M. 
3. 

(3) J (^J (I){t,v)dm2{v)j dmi{t) = J (^J (t>{t,v) dmi{t)^ dm2{v). 

Proof. Using the Lebesgue dominated convergence theorem, it can be 
shown by induction on k that J (f){t, v) dmi{t) exists for v £V and equals 
/ -^(p{t,v) dmi{t). Hence the first assertion. 

Since m2 is of finite order and < Ofc on M x 1/, the function 

t J (j){t,v) dm2{v) is bounded on R. To show it is measurable, consider 
a sequence ipn of C°° functions with compact support that converges, in 
the Schwartz distribution sense, to m2. The sequence ipn can be chosen 
so that there exists a compact K' gV containing the support of ipn for 
all n. Thus, for any t, J (l){t,v)tpniv) dv J (f){t,v) dm2{v) as n — > cxd. Since 
the function (j){t,v)'ipn{'v) is bounded and measurable on M?, the function 
1 1— > / 4>{t^v)ij)n{v) dv is measurable on R for any n. Thus, the function f i— > 
J (j){t,v) dm2{v) is measurable on R as limit of measurable functions. 
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We now show the third assertion. The measure-theoretic version of Fu- 
bini's theorem shows that 

(4) J (^J (j){t,v)iJniv)dv^ dmi{t) = J (^J ^{t,v)dmi{t)^M'")dv. 

By the first assertion, the right-hand side of (4) goes to the right-hand side 
of (3) as n goes to infinity. We prove the third assertion by showing that the 
left-hand side of (4) goes to the left-hand side of (3) as n goes to infinity. 

We say that a set E of C°° functions on M is bounded if the supports 
of all elements of E are contained in the same compact and if there ex- 
ists a sequence Pk such that \C^''\v)\ < (ik for e E, v £R and k>0. It 
is known [20], Theorem V.26, that if a sequence of Schwartz distributions 
converges to a Schwartz distribution, it converges uniformly on any bounded 
set. Let T] he a C°° function on M with a compact support contained in V 
and such that r]{v) = 1 on an open set containing K U K' . Since the set 
E = {(/>(t, .)r], t G M} is a bounded set, / (j){t, v)ipn{v) dv converges uniformly 
to / (j){t,v) dm2{v) as n ^ oo. As a consequence, the left-hand side of (4) 
converges to the left-hand side of (3) as n ^ oo. □ 

Proposition 2.2. ForO<t<T, the function N{^^^^) is C°° with re- 
spect to V. For any k>0 and any closed subset S ofM. — {g(T)}, -^N{^^^^^) 
is bounded for {t, v) £ [0, T) x S. 

Proof. For < t < T, |iiV(^^) = (T - i)-'=/2ArW(H^), where 

A^('^) is the kth derivative of the normal function. Let e > be the minimum 
distance between S and g{T). Since g is continuous at T, there exists s £ 
[0,T) such that \g{t) - g{T)\ < e/2 for s <t <T. For v e S and <t < s, 
'§^^ C'^j'^^l ) is bounded by a constant independent of v and t since the 
function A^^*^) is bounded. On the other hand, for v G S and s < t < T, 
\v - git)\ > e/2. Since the function x''N^''\x) is bounded, \-^N{^^^^)\ < 

( '^^^'^^^^f I A^(^') ( '^^^^^ ) I is also bounded by a constant independent of v and 
t for w G S" and s < t < r. □ 

Proposition 2.3. Ifu<g{0) andv>g{T), then H{u,v) =N{^). 

Proof, li u + Wt > g{T), then t{u) < T a.s. since Wt is a.s. a contin- 
uous function of t. Thus, for v > g{T), 

P{t{u) < T a u + Wt> v)=P{u + Wt>v) 

' u — 



N 



□ 
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Lemma 2.4. 1. For fixed u, the function H{u, ■) is C°° on R — {g{T)}. 
2. If 1/ is a Schwartz distribution with a compact support contained in 
M — {g{T)} then, for any real number u, 

(5) J H{u,v)du{v) = E{U{g{T{u)),T{u);u)). 

Proof. Let mi be the image measure of the random variable t(u) and 
(j){t, v) = N{S^^) if < t < r, and (j){t, v)=0 otherwise. We can write 

H{u, v) = P{t{u) < T a u + Wt>v) 

= ^(l,(„)<rP(n + Wt> v\Tr{u))) 

(6) = E{1,(^^)^tP{Wt -W,i^^)> V - g{T{u))\Triu))) 
= E{(t>{T{u),v)) 

(j){t, v) dmi{t). 

The second equahty follows from the tower law and the fourth one from the 
strong Markov property. It follows from (6), Proposition 2.2 and Lemma 2.1 
that H{u, v) is C°° with respect to v on any open subset 1/ of M such that 
g{T) ^ V. Hence, the function H{u, ■) is C°° on M - {g{T)}. 
We show the second assertion: 

H{u,v)dv{v) = [( I (t>{t,v)dmi{t)] dv{v) 



I 



4){t,v) dv{v) dmiii) 



U{g{t),t;u)dmi{t) 

= E[U{g{T{u)),T{u)-v)). 
The first equality follows from (6) and the second one from (3). □ 

Theorem 2.5. Let fin be a sequence of Schwartz distributions with com- 
pact supports contained in [g(T), +00) such that, for 0<t<T, the sequence 
U{g{t),t; fin) ^ ^ as n —> 00 and, for some constant c, \U{w,t; fin)\ < c for 
all n, <t <T and w < g{t). Then, for u < g{0), U (u, 0; P{t{u) < T) 

as n —> 00 . 

Proof. Given n and e > 0, consider the Schwartz distribution Vn,e such 
that / '4^{v) di'n,e{v) = J ipiv + e) dvn{v) for any C°° function with compact 
support on R. The support of Vn,e is contained is (^(r), +00) and so, by (5), 

(7) / H(u,v)dVn,e{y)=E{U{g{T{u)),T{u);Vn,e)). 
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But it follows from Proposition 2.3 that the left-hand side of (7) is equal to 
U (u, 0; fn,e)- On the other hand, U (x, t; i'n,e) = U{x — e,t; Un)- By letting e go 
to and using the Lebesgue dominated convergence theorem, (7) becomes 

C/(n, 0; Un) = E{U{g{T{u)),T{uy, 

By applying once again the Lebesgue dominated convergence theorem, it 
follows that U {u, 0; Vn) — > P{t{u) < T) as n — > oo. □ 

Corollary 2.6. Let n be a Schwartz distribution with a compact sup- 
port included in [g{T), -\-oo) such that, for some constant c and <t < 
T, U{g{t),t;n) = 1 and \U{w,t; < c if w < g{t) . Then P{t{u) <T) = 
U{u,0;fi) foru<g{0). 

Corollary 2.7. Let n be a a-finite measure on [g{T), +00) such that 
U{g{t),t;fi) = lforO<t<T. Then P{t{u) < T) = U{u,0; fi) for u<g{0). 

Remark 2.8. Let < t < T. A simple time-change argument shows 
that, for u < g{t), the probability that a Brownian motion starting at u 
at time t hits the barrier g between t and T is equal to lim„_>(X) U{u,t; /j-n)- 

3. One-sided last exit time densities. Define 

X = sup{t£{0,T]:Wt = git)} 



and 

Observe that 
Denote by 



a = sup{te{0,T]:Wt>git)}. 
P{a < t) = P{X<tAWT< g{T)). 



(X, ,) = y_ y_ ^-^=== exp ^(l _ ,2) ) du dv 
the cumulative bivariate normal probability distribution function. 

Theorem 3.1. Under the assumptions of Theorem 2.5, for <t <T, 

Furthermore, if g is on an open subset I of [0,T] and the sequence 
U{g{t),t; /j,'^) is uniformly bounded and converges to a continuous function 
of t as t £ L and 00, then 

(9) ±P(a<t) = ±,N'[i^)l^^Um.tM 

for t£l. 
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Proof. Let t e (0,T). By Remark 2.8 

rait) 1 



PiWt < g{t) A A > t)= / ^N' [^] lim U{u, t; du 

/■fW 1 ,fu\ 
= lim / —pN —= ]U(u,t; LLn) du 

rait) 1 



The second equation follows from the Lebesgue dominated convergence the- 
orem and the third one from Lemma 2.1 with 'mi{du) = rn2 = /^n, 

V = Ws.nd<P{u,v) = {l+u')^^N'{^)N{^)l^^g(^t)- Since " 
P{a<t) = P{\< thWt<g{t)) 

= P{Wt < g{t)) - P{Wt < git) A A > t), 
we conclude that (8) holds. On the other hand, since 

'"''(..y;,) = N'(.)N'^-'"^ 



and [19] 



-T^ix,y;p) = -y===N {x)N 



dt \^t' Vt' Vt 

-'(f)(K«l(f 

( 9{t) 



2^t{T-t) VVT^. 
Let J^{t)=jN2{'-§,-^;-^)dfin{v). If [s,t] C / with s < t, 

(11) Ut) - Us) = [ (/ -/I) dO. 

Equation (11) can be derived by applying Lemma 2.1 with mi{d9) = ls<e<t dO, 
m2 = /x„, ^{e,v) = ls<e<t^N2{'-^,-^;-^) and V = R. Combining 
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(10) and (11) yields 

Jn{t) - Jn{s) 



iV'(^) (c/b(«),«;.J^(^) - ^t/b(«),«;.;)) d« 



and so, by the Lebesgue dominated convergence theorem, 
hm Jn{t) - J„(s) 

(12) =n('M-N^'^'^ 



Equation (9) follows by combining (8) and (12). □ 

Example 1. Let fin{v) = 26a+bT + 2be^^(''"''~^'^h{a+bT<v<n} forn > 
and g{t) = a + bt. Then 

as n goes to infinity. By Theorem 2.5, 

(13) P(r(0) < T) = n[^^) + e-^<^^iv(^) 
if a > 0. By Theorem 3.1, for < t < T, 

By symmetry, 

(14) = -L ( '^^^^^^ + 6(iV(6^) - iV(-5VT3i))) 
,//5(i)' 



X iV' 



for < t < T. Note that (13) agrees with the well-known formula [27], Corol- 
lary 7.2.2, for linear barriers crossing. When o = 6 = 0, (14) is identical to 
the well-known first arcsine law [21], page 112. 
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Example 2. /x = ^y^^a = a + &\/r - t. By Corollary 2.7 

(15) p(^r{0)<T) = N^^^/N{b), 

if a + 6\/T > 0. By Theorem 3.1, for < t < T and any values of o, b, 

and 

By symmetry, 
P(A> t) 



and 

^"^^^ - = 2yv(6)yv(-6) Vt(T-t)^ ItT 

for < t < T. Equation (15) agrees with the well-known formula [21], page 
105, for constant barriers crossing when 6 = 0. Here again, (16) is identical 
to the first arcsine law when a = 6 = 0. 

Example 3. /i = \/2Trb5a, where b>T. The conditions of Corollary 2.6 



hold for g{t) = a-^{T -t)ln{^^). Thus P(r(0) < T) = Vb/Tewi-§r) 

L, for any value of a and < t < T, 



a > VTh^{b/T). By Theorem 3.1, for any value of a and < t < T 

9it) 



P{a <t) = N 
and 



Vt 



—P{a<t) = -J— -In- N" ' 



dt ^ - ' 2V t{T-t) T-t \ y/i 



Observe that the function g*{t) = a + y (T — t) ln(7^^) does satisfy the 
equation U{g* {t),t; fi) = 1 but U{w,t;fi) is not bounded for < t < T and 
w < g*{t)- Thus, neither the conditions nor the conclusion of Corollary 2.6 
hold for g*. 
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Example 4. More generally, consider the Hermite polynomial Hn de- 
fined by 

■ exp(-xV2) = (-l)"F„(x) exp(-xV2). 



Thus Hq{x) = 1, Hi{x) = X, H2{x) = — 1 and Hn+i{x) = xHn{x) — H'^{x) 
for n > 0. Moreover, Hn has n distinct zeros, and two successive zeros of Hn 
enclose a single zero of Hn-i- For n > 1, let be the largest zero of Hn and 

let n = V2^b6^J'\ with b > exp(4/2)r"/V^n-i(-Sn)- For < t < T, let x{t) 
be the largest solution to the equation exp(— = (T — t)"^/"^ /b 

and g{t) = a - VT^x{t), with ^(r) = a. By Corollary 2.6, if < ^(0), 
which is equivalent to the conditions 



a>ZnVT and bT~"/^ Hn^ii exp(^-^j < 1, 



,2 



then 

P(r(0) < T) = 6r^"/2^„_i ( ^ 
For any value of a and <t <T, 



expl-- 



dt ^ - ' 2Hn-i{x{t))^t{T -t) V \/t 

When n = 2, (^(t) = a — — tx{t), where x(t) is the largest solution to the 
equation a;exp(— x^/2) = (T — t)/b. The function g{t) is, up to a translation 
and time-reversion, identical to the upper boundary of a symmetric band 
used in [13] as a confidence region for Brownian motion. Whether our tech- 
niques could be used to compute analytically the crossing probability for 
the two-sided band, which has been approximated numerically in [13], is a 
question that deserves further investigation. Note that when n = 2 or n = 3, 
the function g can be expressed in terms of the Lambert function defined 
as the multivalued inverse of the function w i-^ we'^ . This has already been 
observed in [13] when n = 2. Efficient algorithms for computing the Lambert 
function can be found in [3]. 

4. Time-inversion. If Wt is a standard Brownian motion, then 

w - / °' if t = 0, 

~ \tT~^WTyt, if *>0, 

is also ([21], page 21) a standard Brownian motion. By letting g{t) = tT^^g{T'^ /t) 
and 



a = \ni{t>T:Wt>m]. 
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it follows that 

(17) P{3t >T:Wt= m) = Pi^t G (0, T]:Wt= g{t)) 
and, for t>T, 

(18) lp(,<,) = -lp(,<TVt) 
if the latter derivative exists. 

Example 5. It follows from Example 2, (17) and (18) that, for g{t) 
at/T + b^t'^/T-t, 

P{3t >T:Wt = m) = N (^) /Nib) 



if a + h^/T > 0. For t>T and any values of a, b 

d ^ N'(b) ^rifgit) 
—P(& <t) = ^ N' ' ' 



dt ^ - ' 2N{b)ty/t/T -1 V Vt 

5. Two-sided boundaries. Theorems 2.5 and 3.1 can be extended to two- 
sided boundaries as follows. Let go and gi be two continuous functions on 
[0, T] such that gi{t) < go{t) for < t < T. Define 

T{u)=M{te[0,T):u + Wt = go(.t) ov u + Wt = gi{t)} 

and 

a = snp{te{0,T]:Wt>go{t) or Wt<gi{t)}. 

If ^ is a Schwartz distribution or a cj-finite measure, u G M and < t < T, 
let Uo{u,t; fi) = U{u,t; fi) and 



Uiiu,t;f,)= I N(^^=l=jdfi{v). 



Theorem 5.1. Let /io,n [resp. ^i^n) be a sequence of Schwartz distribu- 
tions with compact supports contained in [goiT), +oo) [resp. {—oo, gi{T)]) 
such that, for <t <T, w = go{t) or w = gi{t) , 

1 

^Ui{w,t;fii^n) 1 

i=0 

as n —> oo and, for gi{t) <w< go{t), i € {0, 1}, all n and some constant d , 
\\Ji(w,t;iii^n)\ < c'. Then, for gi{0) <u< goiO), 

1 

P(r(u) < T) = lim ^Ui{u,t; m^n)- 



1=0 
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Furthermore, for <t <T, 
P{(7 < t) 



n— >oo 

i=0 



Finally, if gi is on an open subset I of [0,T], for i € {0, 1}, and the se- 
quence U{gi{t),t; fiQ ^ — fj^'i^) is uniformly bounded and converges to a con- 
tinuous function of t as t £ I and n — > oo, then 

for t€l. 

Example 6. Let go{t) = a and gi{t) = b for < t <T, where b < a. 
Define fj,o,n = '2J2]=o{-iy 6a+j{a-b) and ;Ui,„ = 2X;"=o(-iy^b+j(6-a)- Then 

t/.(..,.,„) = 2t(-iyAr( '-^»'"-^+^-(^-°' ) 

is positive and upper-bounded by 1 for i G {0, 1} and b <w < a. If w = a or 
w = b, then 

i:f/.(».t;,.„J = 1 + 2(-l)..A,((!i±|£^) ^ 1 

as n ^ oo. If 6 < < a, Theorem 5.1 yields the well-known formula (see, 
e.g., [21], Exercise 3.15, page 111, for a direct derivation) 

P(r(0)<T)=2f:(-lFiv( -° + ^;<^-°> ) 

Furthermore, for < t < T, 

dt ^ - ' ^2TTt{T -t)\ \VtJ \Vt 



l-h2^(-iyexp 
j=i 



2{T-t) 
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Example 7. Let /xq = c ^Sa and /ii = c ^5b, where b < a and 2N{^^) < 
c < 1. For < t < T, let go{t) be the largest solution to the equation 

and gi{t) = a + b — go{t). If 6 < < a, then, by Theorem 5.1, 

PWO)<r,^c-.iv(^)..-.;v(-^) 

and 

6. Schwartz distributions and the method of images. Let Bt, t>0, be 
a Brownian motion and / a continuous function on an interval [0,T'] such 
that /(O) > 0. Define the first hitting time r by 

T = mf{te[0,T'):Bt = f{t)}. 

Fix a real number T £ (0, T'] and let h be the function from M to M"*" defined 
by 

h{u)=P{T<T A Bt<u). 

For < f < T, let g{t) = f{T - 1) and Wt = Br^t - Bt- Thus, Wt is a Brow- 
nian motion [21], page 21, on [0,T] and = —Bt- As before, the random 
variable t{u) and the function H are defined via (1) and (2), respectively. 
The following proposition shows formally that the probability that a Brow- 
nian motion starting at u at time and ending at at time T hits g equals 
the probability that a Brownian motion starting at at time and ending 
at u at time T hits /. 

Proposition 6.1. Foru< g{0), h'{u) = -^{u,0). 

Proof. By Lemma 2.4, H{u, •) has a first derivative with respect 
to u at V = since g{T) = /(O) > 0. Similarly, using the relation h{u) = 
E{lr<TN{^^^^)), it can be shown that /i is C°° on M - {^(0)}. For e>0, 

P{t{u) <TA 0<u + WT<e) 
(19) = P(r(n) <rA u-£<Bt<u) 

>P{t<TA u-£<Bt<u). 
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To show the second equation, we observe that if r < T and u — e < Bt < u, 
then Bt = fit) for some t G (0,r). Thus, Wr^t + Bt = g{T - t) and so 
Wt-i + u> g{T — t). Since u < g{0) and W is a.s. continuous on [0,T], we 
conclude that, a.s., Wf + u = g{t') for some t' £ [T — t, T). Hence t{u) < T 
a.s. It follows from (19) that H{u, 0) — H{u, e) > h{u) — h{u — e). By letting 
e ^ 0, we conclude that h'{u) < —^{u,0). 
Similarly, 

P{t{u) < T a -e<u + WT<0) 

(20) = P(r(n) <rA u<BT<u + e) 

<P(r<TA u< Bt <u + e). 

We show the second equation by noting that if t{u) < T and u < Bt <u + e, 
then u + Wt = g{t) for some t G (0, T). Thus, u + Bt^i - Bt = f{T - t) and 
Bt-i > f(T — t). Since /(O) > 0, we conclude that, a.s., Bf = f(t') for some 
t' G (0, r - t], and so r < T. It follows from (20) that H{u, -e) - H{u, 0) < 
h{u + e) - h{u), and so h'{u) > -|^(n,0). Thus h'{u) = -|^(u,0), as de- 
sired. □ 

Theorem 6.2. Let ^ he a Schwartz distribution with compact support 

contained in the interval (/(0),+oo). Assume that / exp( ^-^^^2f~" )^/^('^) = 1 
forte (0,r']. Then 

(21) P{r < T) = N^=^^ + U{f{T)A /.). 
If f is on an open subset I of [0, T'], then, for T £ I, 

Proof. Consider the Schwartz distribution v = 6q — n' . For < t < T, 
Uig{t),t;u) = U{f{T -t),t;5'o- fi') 

= 0. 

It follows from Lemma 2.4 that j H{u,v)du{v) = for any real number u 
and so, by Proposition 2.3, 




= -C/(tx,0;/i'). 
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Using Proposition 6.1, we conclude that h'{u) = U{u,0;ij,') for u < f{T). 
Since h is a continuous function of u and is on {—oo,f{T)), it follows 
that, for u < f{T), 



(24) 



h{u) = I h'{x) dx 



= U{u,0;fi). 

The third equation follows from Lemma 2.1 with mi{dx) = j^^, ^-2 = 
(j){x,v) = {l + x^)N'{^)l^^u andy = (/(0),+oo). 

We infer from (24) that P{t<T ABt< f{T)) = U{f{T),0; n). On the 
other hand, since the events r < T A Bt > f{T) and Bt > f{T) coincide up 
to a negligible set, P{t <T ABt> f{T)) = N{-f{T)/VT). This concludes 
the proof of (21). 

Assume now that / is on an open subset / of [0,r'] containing T. 
Then 

d {fiT)\,r>ffiT)\ 



27^3/2 

The first equality can be justified by fixing a real number s such that [s, T] C 
/ and applying Lemma 2.1 with mi{dt) = ls<t<Tdt, m2 = 1^ = M and 

^{t,v) = ls<t<TiNii^). □ 



-2 



Example 8. ^ = bd'^, where a > and b > aexp(l — ^). Let /(O) = a/2 
and f(t) = a + ty[t)/a for < t < T', where y{t) is the smallest solution to 

2 

the equation y{t) exp{y{t)) = — aexp(— |j)/6. The function / can be easily 
expressed in terms of the Lambert function. By Theorem 6.2, it follows after 
some calculations that, for <T <T' , 
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and 

7. Concluding remarks. We conclude with tiie following remarks: 

1. Our results can be adapted to other Markov processes, including discrete- 
time Markov chains. 

2. The condition \U{w,t; < c for all n, < t < T and w < g{t) in The- 
orem 2.5 can be relaxed as follows: \U{w,t] < c for all n, s <t <T 
and w < g{t), for some fixed s G [0,T). 

3. Another approach to prove Theorem 2.5 is to show that U{u, t; ^) satisfies 
the partial differential equation 

dU Id'^U _ 

and to use the martingale properties of U{Wt,t; fj.). 

4. The Schwartz distributions in our examples were essentially found by in- 
spection. The questions of establishing necessary and sufficient conditions 
for the existence of a sequence of Schwartz distributions satisfying the 
conditions of Theorem 2.5 for a given function g, and determining such 
a sequence if it exists, are left for future research. A numerical method 
for calculating cr-finite measures associated to a given boundary via the 
classical method of images has been used in [6, 16] to approximate the 
first hitting time density in cases where the exact density is not known 
or difficult to compute. 

Acknowledgments. The author would like to thank Monique Jeanblanc, 
Eric Sere and an anonymous referee for constructive comments. 
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